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I Abstract. We review recent results on the inclusive electroproduction of light hadrons at next-to- 

. leading order in the parton model of quantum chromodynamics implemented with fragmentation 
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^ ■ 1. INTRODUCTION 

O ■ 

^ I In the framework of the parton model of quantum chromodynamics (QCD), the inclusive 

^ I production of single hadrons is described by means of fragmentation functions (FFs) 

lO ■ ^ai^il^)- At lowest order (LO), the value of /i) corresponds to the probability 

^ ! for the parton a produced at short distance l//i to form a jet that includes the hadron 

^ I h carrying the fraction x of the longitudinal momentum of a. Analogously, incoming 

hadrons and resolved photons are represented by (non-perturbative) parton density func- 
Qh! tions (PDFs) F^//,(x,jU). Unfortunately, it is not yet possible to calculate the FFs from 

^ I first principles, in particular for hadrons with masses smaller than or comparable to the 

■ asymptotic scale parameter A. However, given their x dependence at some energy scale 

• T-j , jU, the evolution with /i may be computed perturbatively in QCD using the timelike 

I Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. Moreover, the factor- 

ed ' ization theorem guarantees that the D^(;c, /i) functions are independent of the process 

in which they have been determined and represent a universal property of h. This enti- 
tles us to transfer information on how a hadronizes to in a well-defined quantitative 
way from e+e^ annihilation, where the measurements are usually most precise, to other 
kinds of experiments, such as photo-, lepto-, and hadroproduction. Recently, FFs for 
light charged hadrons with complete quark flavour separation were determined through 
a global fit to e+e data from LEP, PEP, and SLC [1] thereby improving a previous 
analysis [2]. 

The QCD-improved parton model should be particularly well applicable to the in- 
clusive production of light hadrons carrying large transverse momenta (pr) in deep- 
inelastic lepton-hadron scattering (DIS) with large photon virtuality (Q^) due to the pres- 
ence of two hard mass scales, with Q^.Pj ^ A^. In Fig. 1, this process is represented 
in the parton-model picture. The hard- scattering (HS) cross sections, which include col- 
ored quarks and/or gluons in the initial and final states, are computed in perturbative 
QCD. They were evaluated at LO more than 25 years ago [3]. Recently, the next-to- 




leading-order (NLO) analysis was performed independently by three groups [4, 5, 6]. 
A comparison between Refs. [5, 6] using identical input yielded agreement within the 
numerical accuracy. 

The cross section of e^p e^K^ +X in DIS was measured in various distributions 
with high precision by the HI Collaboration at HERA in the forward region, close to 
the proton remnant [7, 8]. This measurement reaches down to rather low values of 
Bjorken's variable xb — / {2P-q), where P and q are the proton and virtual-photon 
four-momenta, respectively, and = —q^, so that the validity of the DGLAP evolution 
might be challenged by Balitsky-Fadin-Kuraev-Lipatov (BFKL) dynamics. 

In Ref. [5], the HI data [7, 8] were compared with NLO predictions evaluated with 
the KKP FFs [2]. In Section 2, we present an update of this comparison based on the 
new AKK FFs [1]. Our conclusions are summarized in Section 3. 



2. COMPARISON WITH HI DATA 

We work in the modified minimal-subtraction (MS) renormalization and factorization 
scheme with «/ = 5 massless quark flavors and identify the renormalization and factor- 
ization scales by choosing jl^ = ^ [Q^ -f- (pj)^]/2, where the asterisk labels quantities in 
the 7*p center-of-mass (cm.) frame and is varied between 1/2 and 2 about the default 
value 1 to estimate the theoretical uncertainty. At NLO (LO), we employ set CTEQ6M 
(CTEQ6L1) of proton PDFs [9], the NLO (LO) set of AKK FFs [1], and the two-loop 

(one-loop) formula for the strong-coupling constant ai"^^(ju) with A^^) = 226 MeV 
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FIGURE 2. HI data on (a) da/dp*j, (b) da/dxE, and (c) do/dxB for 2<Q^ < 4.5 GeV^, 4.5 <Q^ < 
15 GeV^, or \5<Q^ < 70 GeV^, and on (d) do/dQ^ from Refs. [7] (open cii-cles) and [8] (solid cii-cles) 
are compared with our default LO (dashed histograms) and NLO (solid histograms) predictions including 
theoretical uncertainties (shaded bands). The QCD-correction (K) factors are also shown. 



(165 MeV) [9]. 

The HI data [7, 8] were taken in DIS of positrons with energy = 27.6 GeV on 
protons with energy Ep = 820 GeV in the laboratory frame, yielding a cm. energy of 
^/S = l^EeE'p = 301 GeV The DIS phase space was restricted to 0.1 < y < 0.6 and 
2 < 2^ < 70 GeV2, where y = Q^/{xbS). The TZ*^ mesons were detected within the 
acceptance cuts pj > 2.5 GeV, 5° < < 25°, and xe > 0.01, where 6 is their angle 
with respect to the proton flight direction and E = xsEp is their energy in the laboratory 
frame. The comparisons with the our updated LO and NLO predictions are displayed in 
Figs. 2(a)-(d). 



3. CONCLUSIONS 



We calculated the cross section of ep eir +X in DIS for finite values of pj at LO 
and NLO in the parton model of QCD [5] using the new AKK FFs [1] and compared it 
with a precise measurement by the HI Collaboration at HERA [7, 8]. 

We found that our LO predictions always significantly fell short of the HI data 
and often exhibited deviating shapes. However, the situation dramatically improved 
as we proceeded to NLO, where our default predictions, endowed with theoretical 
uncertainties estimated by moderate unphysical-scale variations, led to a satisfactory 
description of the HI data in the preponderant part of the accessed phase space. In other 
words, we encountered K factors much in excess of unity, except towards the regime of 
asymptotic freedom characterized by large values of pj and/or Q^. This was unavoidably 
accompanied by considerable theoretical uncertainties. Both features suggest that a 
reliable interpretation of the HI data within the QCD-improved parton model ultimately 
necessitates a full next-to-next-to-leading-order analysis, which is presently out of reach, 
however. For the time being, we conclude that the successful comparison of the HI 
data with our NLO predictions provides a useful test of the universality and the scaling 
violations of the FFs, which are guaranteed by the factorization theorem and are ruled 
by the DGLAP evolution equations, respectively. 

Significant deviations between the HI data and our NLO predictions only occurred 
in certain comers of phase space, namely in the photoproduction limit 0, where 

resolved virtual photons are expected to contribute, and in the limit 77 — > oo of the 
pseudorapidity 77 = — ln[tan(0/2)], where fracture functions are supposed to enter the 
stage. Both refinements were not included in our analysis. Interestingly, distinctive 
deviations could not be observed towards the lowest xg values probed, which indicates 
that the realm of BFKL dynamics has not actually been accessed yet. 
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